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ABSTRACT 

The  SANBAR  barotropic  hurricane  forecast  track  model  is 
applied  to  the  four  named  Atlantic  tropical  vortices  from  8 
September  to  13  September  1971.   All  segments  of  the  program, 
one  of  the  three  major  objective  tools  used  at  the  National 
Hurricane  Center  during  1971,  are  described.   The  data  base 
includes  NHC's  history  tape  plus  subjective  use  of  film  strips 
from  the  ATS  III  satellite. 

SANBAR  forecasts  are  plotted  and  compared  with  subsequent 
storm  or  hurricane  movements.   The  effects  of  storm  parameter 
variation  are  discussed  along  with  an  NHC  modification  to 
SANBAR  promising  improved  results. 
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INTRODUCTION 

Since  the  decade  of  the  fifties,  attempts  have  been  made  to 
provide  hurricane  forecasters  with  a  viable  objective  hurricane 
track  forecast.   During  the  1971  hurricane  season,  the  National 
Hurricane  Center  (NHC)  had  three  primary  objective  forecast  tools 
from  which  to  draw  before  releasing  its  official  advisories  on 
tropical  storm  and  hurricane  movements.   The  first,  called  NHC-67, 
is  a  statistical  method  developed  by  Killer  et  al  (1968)  which 
uses  a  multiple  screening  regression  technique  dependent,  now, 
upon  prognostic  and  analyzed  values  of  pressure-height  data  as  well 
as  past  storm  motion  and  other  synoptic  parameters.   The  second, 
known  as  HURRAN  (for  hurricane  analog)  produced  by  Hope  and  Neumann 
(1970)  selects  tracks  of  similar  hurricanes  dating  back  as  far  as 
1886  for  possible  application  to  the  existing  tropical  storm  or 
hurricane.   The  third  method,  a  dynamic  one  called  SAN3AR  (for 
Sanders  barotropic) ,  developed  by  Professor  Frederick  Sanders  of 
the  Massachusetts  Institute  of  Technology  is  the  subject  of  this  work, 
This  model  is  based  upon  the  hypothesis  that  the  advection  of  mean 
tropospheric  vorticity  in  the  air  column  containing  the  storm  pri- 
marily governs  the  motion  of  the  cyclone.   The  depth  of  the  appro- 
priate air  column  is  taken  to  be  from  1000  to  100  mb. 

Following  attempts  using  vortex  separation  (Kasahara,  1957),  the 
technique  of  barotropic  calculation  without  separation  of  the  storm 


vortex  from  the  basic  flow  has  been  investigated  by  Sanders  (1959) , 
Birchfield  (1960,  1961),  King  (1966),  Ahn  (1967),  as  well  as 
Sanders  and  Burpee  (1968).   Next  came  the  substitution  of  the  pre- 
vious manual  analysis  with  an  automated  analysis  scheme  (Sanders, 
1968,  1970).   The  current  report  investigates  the  results  of  SANBAR 
applied  to  the  several  cyclones  existing  in  the  tropical  Atlantic 
during  the  period  September  8-13,  1971.   An  effort  is  made  to  im- 
prove the  wind  data  available  from  the  NHC  history  tape,  and  use 
is  made  of  satellite  information  in  regions  where  soundings  are  un- 
available.  Additionally,  storm  parameters  are  varied  to  display 
the  effects  of  such  variations  upon  the  forecast  tracks. 
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DATA  PREPARATION 

Since  previous  investigations  of  King  (1966)  and  Ahn  (1967) 
indicate  that  the  10  mandatory  levels  represent  an  optimum  vertical 
sample  from  which  to  work,  the  average  wind  for  each  reporting 
station  is  approximated  by  the  trapezoidal  rule: 


V  -£t£<*£i«0j/i  +  /4-y^+jfcfcY, 


where:    p  =  lOOOmb,  p,  =  850mb,  .  .  ,pn  =  lOOmb 
o  1  9 

V.=  the  horizontal  wind  vector  for 
the  level  specified 
(This  and  all  subsequent  computations  were  performed  on  the  IBM 
370  Model  155  Computer  at  the  Information  Processing  Center  of 
M.I.T.) 

Winds  for  the  10  mandatory  levels  were  obtained  from  magnetic 
history  tapes  made  available  by  NHC.   For  the  sample  period  chosen 
(O000Z  September  6,  1971  through  1200Z  September  13,  1971)  the 
number  of  soundings  available  at  a  specific  synoptic  time  ranged  from 
126  to  146  within  the  area  of  the  grid  —  to  be  described  later.   To 
improve  the  data  base  so  supplied  by  the  tape,  teletype  reports 
were  scanned,  and  missing  soundings  added  when  found.   When  interior 
level  winds  were  missing,  they  were  interpolated  from  surrounding 
levels;  winds  missing  from  the  bottom  or  the  top  were  computed  by 
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extrapolation  of  the  lowest  or  highest  level  wind  available,   Finally, 
those  soundings  having  fewer  than  two  reported  winds  from  the  lower 
A  levels  or  fewer  than  two  from  the  upper  6  levels  were  rejected 
and  thus  considered  to  be  missing.   Using  this  convention  on  the 
16  synoptic  times  involved,  the  number  of  acceptable  soundings  then 
ranged  from  100  to  122.   It  may  be  of  interest  to  note  that  the  number 
of  soundings  usable  south  of  30°N  was  limited,  having  a  median  of 
24  while  ranging  from  22  to  27. 
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OBJECTIVE  ANALYSIS 

A.   Development  Program 

After  the  observed  mean  winds  are  prepared,  the  SANBAR  forecast 
entails  the  execution  of  3  computer  programs.   The  first,  called 
the  "development"  program  may  be  executed  far  in  advance  of  the  2 
operational  programs. 

The  development  program,  written  by  Mr.  Peter  P.  Chase  of  the 
National  Hurricane  Research  Laboratory  following  the  procedures  pro- 
posed by  Eddy  (1967),  provides  regression  equations  which  will  be 
used  to  compute  winds  for  missing  stations  as  well  as  winds  for  all 
2655  grid  points.   The  numerical  grid  consists  of  a  45  x  59  array 
of  grid  points  extending  from  the  equator  to  55°N  Latitude  and  from 
36.5°W  to  123. 5°W  Longitude  on  a  Mercator  projection  true  at  22.5°N. 
The  mesh  length  is  154  km  (at  22.5°N).   Fundamentally,  this  program 
is  based  on  an  estimate  of  the  correlation  of  the  observations,  or 
more  specifically  the  observed  departures  of  zonal  and  meridional 
wind  components  from  their  respective  latitudinal  means,  as  a  function 
of  distance  between  observing  points.   Figure  1  shows  this  estimated 
correlation  based  upon  the  1713  input  soundings.   However,  stations 
within  the  influence  radius  of  cyclones  are  excluded  from  the  com- 
putation, as  are  stations  having  wind  greater  than  5  standard  devia- 
tions from  the  mean  of  the  specific  synoptic  time  period  or  from  the 
latitudinal  mean.   (Stations  having  wind  greater  than  3  standard 
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Fig.  1  —  Correlation  as  a  function  of  separation 
distance  for  departures  of  vertically  averaged  wind 
from  zonal  average  values 


deviations  from  these  means  are  "flagged"  to  permit  investigation  by 
operating  personnel.)   Contrasted  with  a  previous  sample  of  soundings 
which  computed  the  correlation  dropping  to  zero  at  700-800  nautical 
miles  (Sanders,  1970),  in  the  current  sample  the  correlation  drops  to 
zero  at  560-580  nm.   This  would  seem  to  imply  that  within  the  one 
week  period  of  this  sample,  the  prevailing  size  of  the  circulation 
was  smaller  than  that  of  the  previous  study  —  the  previous  sample 
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in  particular  did  consist  of  data  from  several  hurricane  seasons. 
The  current  sample  also  displays  higher  correlation  between  stations 
less  than  175  nm  apart  and  lower  correlation  for  stations  greater 
than  175  nm  apart  than  did  the  previous  sample. 

Table  1  presents  a  sample  of  the  regression  equations  generated 
by  the  development  program  —  equations  that  are  used  in  the  following 
"bperations"  program  should  a  station  be  missing  during  a  specific 
synoptic  time.   Several  stations  existing  outside  the  grid  are  included 
with  the  data  since  they  may  be  used  in  the  regression  equations  (i.e., 
they  are  within  580  nm) .   Similar  regression  equations  are  generated 
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Table  1.   Regression  equation  coefficients  for  computing 
winds  for  missing  stations.   (Winds  are  computed  for  the 
station  in  the  left-most  column  using  wind  components  of 
stations  to  the  right  with  the  associated  coefficients  — 
zonal  component  on  top  row,  meridional  component  on  bottom 
row.   Unexplained  variance  for  both  the  u  and  v  components 
are  shown  also.) 
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to  provide  winds  at  each  grid  point  —  winds  which  are  subsequently 
used  in  forecast  calculations. 

Required  as  input  to  the  development  program  are: 

1.  A  station  list  supplying  block  and  index  numbers  as  well  as 
latitude  and  longitude  (the  185  northern  hemisphere  stations  in  this 
investigation  included  several  now  defunct  stations  should  SANBAR  be 
tested  on  storms  from  recent  previous  years); 

2.  Mean  winds  from  available  soundings  (in  this  study,  1713 
soundings  for  the  16  synoptic  time  periods  during  the  sample  period 
chosen) ; 

3.  The  tropical  storm  and  hurricane  data  for  each  synoptic  time 
period  (latitude,  longitude,  and  maximum  influence  distance). 

In  deriving  the  correlation  curve  and  regression  equations,  the 
statistics  of  the  wind  data  are  assumed  to  be  isotropic  (aside  from 
the  latitudinal  trend),  homogeneous,  and  stationary.   As  noted  by 
Sanders  (1968) ,  the  failure  of  the  data  to  meet  these  requirements 
does  not  appear  to  lead  to  crucial  errors  in  the  final  analysis.   We 
have  noted  in  particular  an  instance  of  non-stationarity  in  the  com- 
parison of  samples  discussed  above.   Nevertheless,  this  is  not  meant 
to  imply  that  continued  research  in  this  area  is  not  welcomed. 

i 
B.   Operations  Program 

! 

The  operations  program,  also  written  by  Mr.  Chase,  must  be  executed 
each  time  in  preparation  for  the  barotropic  forecast.   Its  ultimate 
purpose  is  to  provide  u  and  v  wind  components  for  each  grid  point  in 
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the  array.   Required  as  input  to  the  program  are: 

1.  Mean  winds  for  soundings  at  the  specific  synoptic 
time  for  which  a  forecast  is  desired  (  Additional  esti- 
mated mean  winds  must  be  provided  for  data-sparse  regions, 
as  described  below.); 

2.  Hurricane  parameters  for  the  specific  synoptic  time 
(latitude,  longitude,  maximum  influence  distance,  maximum 
surface  wind,  and  eye  diameter  —  the  maximum  influence 
distance  is  subjectively  estimated  from  the  surface  analysis). 

Mean  winds  from  soundings  are  supplemented  by  "bogus"  winds  for 
the  44  locations  indicated  over  data  sparse  and  oceanic  regions  shown 
in  Figure  2.   Below  35°N,  the  first  estimate  of  the  mean  wind  at  bogus 
points  was  made  from  200  mb  and  surface  analyses  prepared  by  NHC. 
Assuming  a  linear  variation  of  the  wind  from  the  surface  to  100  mb, 
the  contributions  of  the  200  mb  and  of  the  surface  winds  are  .56  and 
.44  respectively.  However,  comparison  of  mean  winds  already  computed 
during  the  data  preparation  for  certain  oceanic  stations,  such  as 
Bermuda,  often  indicated  that  weighting  factors  should  be  near  .30 
and  .70  respectively  —  i.e.,  a  heavier  weight  given  to  the  surface 
wind  than  to  the  200  mb  wind.   Accordingly  the  weight  given  the  200  mb 
wind  and  the  surface  wind  was  subjectively  determined.   A  great  assist 
in  preparing  the  bogus  winds  for  0O00Z,  came  from  viewing  film  strips 
from  the  Applications  Technology  Satellite  III  (ATS  III). 
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Fig.  2  —  SANBAR  Numerical  Grid  (X's  indicate 
the  forty-four  bogus  points.) 


Above  35°N  it  is  found  that  the  geostrophic  wind  measured  from  the 
500  mb  analysis  provides  a  very  close  estimate  of  the  mean  wind. 

Immediately  after  reading  in  the  -station  and  bogus  data,  the 
input  parameters  for  any  storms  present  are  used  to  identify  those 
winds  which  are  within  the  maximum  influence  radius  of  the  storm. 
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Within  this  radius,  the  storm/hurricane  circulation  is  subtracted  from 
the  observed  mean  winds  leaving  a  "recalculated"  or  storm-purged  wind. 
This  model  hurricane  wind  determined  empirically  as 

£q0£_        IS 

where:  V    =  Maximum  surface  wind 
max 

R    =  Distance  from  storm  center 

R    =  Max.  influence  distance 
max 

E    =  Hurricane  eye  diameter 
is  vectorially  subtracted  from  the  input  wind  in  an  attempt  to  depict 
what  the  undisturbed  wind  would  have  been.   Then  anomalous  winds  are 
rejected  as  in  the  development  program.   If  winds  from  stations  used 
in  grid-point  regression   equations  are  missing,  they  are  reconstructed 
using  station  regression  equations  from  the  development  program.   How- 
ever, should  this  be  impossible  because  of  additional  missing  data, 
ad  hoc  regression  equations  are  generated  from  the  winds  which  are 
present. 

The  analysis  for  the  entire  array  is  then  carried  out  by  solving, 
at  each  grid  point,  the  regression  equation  for  the  zonal  and  meridi- 
onal components  of  the  wind.   Finally  at  those  grid  points  within 
the  influence  radius  of  a  storm,  the  zonal  and  meridional  components 
of  the  storm  are  added  to  the  values  estimated  from  the  regression 
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equations  to  complete  the  analysis. 

Additional  computations  are  made  to  show  the  closeness  of  fit  of 
the  analysis  to  the  data,  but  they  will  not  be  discussed  here. 
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FORECAST  MODEL 

The  Sanders  barotropic  hurricane  forecast  track  program,  written 
with  the  programming  assistance  of  Dr.  Robert  W.  Burpee  and  Mr. 
Frederick  S.  Zbar  while  they  were  graduate  students  at  M.I.T.,  has 
at  its  beginning  a  brief  section  entitled  Diagnosis.   Using  the  opera- 
tions program  wind  analysis  components  as  input,  relative  vorticity 
is  computed  by  a  simple  centered  finite-difference  approximation 

where:  %  represents  relative  vorticity 

jU_  represents  the  eastward  wind  component 

/IT'  represents  the  northward  wind  component 

/±  represents  the  mesh  length 

JT  represents  the  earth's  radius 

Q  represents  latitude 

sjy^  represents  the  map  factor 

J.  J,         represent   the  x  and  y  grid  positions 
respectively 

for  all  interior  grid  points,  with  boundary  values  set  to  zero.   After 

smoothing  the  relative  vorticity  field,  stream  function  values  are 

obtained  by  solving  the  Poisson  equation  using  a  9-point  laplacian 

operator  and  an  extrapolated  Liebmann  relaxation  technique 


T*4        V       L  v    *    W    "V'    V'' 

■f  (A  .  *#  .-f^.  .-f^    -  ga*-  / 
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44 


where:  t/   represents  the  iteration  step 
W     represents  the  stream  function 

J 

•^      represents  smoothed  rel.  vorticity 
C.  and  C     represent  coefficients  including 
overrelaxation 
with  the  Neumann  boundary  condition.   Thus  the  initial  values  of  rela- 
tive vorticity  and  stream  function  are  prepared  for  use  by  the  forecast 
program. 

The  forecast  program  begins  by  calculating  the  absolute  vorticity 


%  =  </  +  *i 


where:  T\       represents  the  absolute  vorticity 
«f        represents  the  Coriolis  parameter 
over  the  entire  grid.   After  the  assumption  of  vorticity  conservation, 
the  finite-difference  Jacobian  operator  is  used  to  obtain  an  approxi- 
mation of  the  horizontal  advection  of  absolute  vorticity.   Then  follows 
a  solution  for  the  change  of  stream  function  with  time  by  the  itera- 
tion formula 
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where:  C„  and  C,  are  coefficients  including  overrelaxation 

V   represents  the  horizontal  velocity  of  the 
mean  wind 

^   represents  the  horizontal  del  operator 

using  the  Liebmann  sequential  relaxation  technique,  with  the  outer  2 

rows  held  constant.   After  the  first  step,  normally  4-9  iterations 

are  required  for  the  solution  of  each  succeeding  time  step.   Conver- 

2 
gence  is  accepted  when  an  absolute  value  change  of  less  than  1  m 

-2  .    L   .   , 
sec   is  obtained: 
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Then  using  a  centered  time  step  (except  on  the  first  step,  of  course) 
the  succeeding  value  of  the  stream  function  is  computed: 

(P.    .  /    -  (b         L     +    2.   At(>L      \ 


where:   k  represents  the  number  of  the  time  stop 
£b    represents  the  time  increment  (30  min.) 
Thereafter,  the  successive  values  of  absolute  vorticity  are 
computed  by  the  finite-difference  formula: 
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[Readers  desiring  a  reference  text  for  applications  of  finite  dif- 
ference techniques  to  meteorology  are  referred  to  the  recent  publi- 
cation by  Haltiner  (1971).] 

At  the  end  of  time  steps  24,  48,  72,  96,  120  and  144  (i.e., 
each  twelve  hours),  a  smoothing  subroutine  is  applied  to  the  field 
of  stream  function  —  and  a  smoothing  of  the  vorticity  field  is 
thereby  indirectly  accomplished. 

The  tropical  storm  is  identified  at  each  printout  time  by 
reference  to  the  associated  stream  function  minimum  and  vorticity 
maximum.   As  the  two  drift  apart  during  the  forecast,  due  to  the 
effects  of  truncation  error,  the  convention  has  been  established  to 
position  the  storm  halfway  between  them.   At  times,  the  minimum  in 
the  stream  function  disappears  during  the  course  of  the  forecast  — 
occurring  often  when  the  storm  is  weak  or  superimposed  on  an  ever- 
stronger  basic  current.   In  this  event,  the  position  of  the  storm  is 
presumed  to  be  indicated  by  the  vorticity  maximum  alone. 
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ATLANTIC  HURRICANES  AND  TROPICAL  STORMS  (Sept.  8-13,  1971) 

The  period  September  8-13,  1971  presented  a  most  profitable  oppor- 
tunity for  application  of  the  SANBAR  model.   Not  only  were  storms  loca- 
ted in  the  Caribbean  Sea,  Gulf  of  Mexico,  and  the  western  Atlantic, 
but  never  were  there  fewer  than  two  named  storms.   Hurricane  Nanette 
located  off  Baja  California  in  the  Pacific  Ocean  was  included  in  the 
analysis  from  0000Z  September  8th  through  1200Z  September  9th  to  add 
to  the  realism  of  the  mean  flow  pattern. 

In  preparation  for  studying  the  forecasts,  it  should  be  noted 
that  the  tracks  plotted  for  the  actual  storms  are  those  published 
after-the-fact  by  NHC  (Simpson  and  Hope,  1972).   However,  the  forecasts 
start  from  the  real-time  operational  positions.   Differences  in  the 
actual  track  positions  and  forecast  initial  positions  often  are 
further  increased  by  the  mid-point  between  the  vorticity  maximum 
and  the  stream  function  minimum  being  removed  from  the  operational 
position  at  time  zero. 

Resisting  temptation  —  and  thus  saving  the  best  for  last  —  I 
will  present  the  storms  in  chronological  order. 
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EDITH 

In  the  forecasts  for  Edith  (Figure  3)  the  first  two  forecasts 
have  creditable  accuracy  for  landfall  near  the  Honduras/Nicaragua 
border,  but  curvature  toward  the  north  is  not  hinted  until  the 
090000Z  forecast.   (Table  2  gives  the  relevant  storm  parameters 
used  for  Edith's  forecasts.) 


Table  2.   Storm  Parameters  for  Hurricane  Edith   (Eye 
diameter  =  20  nm.  throughout) 


Forecast  time 

080000Z 
081200Z 
090000Z 
0912OOZ 
100000Z 
101200Z 
110000Z 
111200Z 
120000Z 
121200Z 
130000Z 


Max.  sfc.  wind  (knots) 

61 

66 

92 
132 
122 

66 

70 

52 

35 

35 

35 


Influence  radius (nnp 

350 

300 

350 

350 

200* 

350 

400 

325 

300* 

300* 

300 


*  indicates  re-run 

By  the  091200Z  forecast,  Edith  was  embedded  in  a  northwestward  mean 
current,  on  the  western  periphery  of  the  subtropical  ridge  extending 
westward  over  the  Caribbean.   The  direction  of  the  forecast  was  ex- 
cellent, but  the  forecast  speed  was  slow  as  we  will  come  to  expect. 
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The  100000Z  forecast  as  shown  is  a  re-run.   The  original  forecast 
is  not  shown,  because  any  operator,  having  noticed  that  its  12-hour 
forecast  reversed  Edith's  course  into  the  strong  easterlies,  could 
have  taken  corrective  action.   In  fact,  an  overly  zealous  hurricane 
influence  radius  (400  nautical  miles)  had  been  used  as  a  storm  para- 
meter on  the  original  operational  analysis  program.   Since  Edith 
was  at  the  time  a  dangerously  strong  hurricane  (122  knots  maximum 
surface  wind),  the  subtracting  of  the  model  hurricane  wind  from  the 
mean  wind  reported  at  Swan  Island,  station  78501,  (097°/27.6  knots) 
left  a  strong  westerly  wind  at  the  station,  (294/15.9).   Accordingly, 
a  quick  re-run  of  the  analysis  program  using  a  200  nm.  radius  for 
the  hurricane  —  thereby  greatly  reducing  the  model  storm  wind  at 
Swan  Island  —  provided  a  new  recalculated  (storm  purged)  wind  at 
Swan  Island  of  090/13.1.   Therefore,  the  many  grid-points  nearby 
which  use  the  Swan  Island  wind  in  their  regression  equations  were 
given  more  believable  easterly  winds  before  the  addition  of  the 
model  hurricane  wind.   (Again,  although  the  forecast  speed  was  slow, 
the  72-hour  error  was  only  140  nm.  —  comparable  with  the  best 
state-of-the-art  24-hour  error.) 

The  101200Z  forecast  was  poor,  caused  by  a  building  westward 
of  the  subtropical  ridge  by  the  baro tropic  forecast;  however,  the 
110000Z  recovered  to  have  only  a  60  nm.  error  at  72  hours,  despite 
its  poor  initial  direction. 

Commencing  at  111200Z  the  deterioration  of  SANBAR's  accuracy 
in  forecasting  for  Edith  becomes  obvious.   Although  the  first  12-hour 
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forecast  is  quite  good,  thereafter  Edith  is  taken  erroneously  north- 
eastward around  the  subtropical  high. 

The  presence  of  tropical  storm  Fern  (to  be  discussed  next),  going 
inland  over  Corpus  Christi,  Texas,  may  have  been  instrumental  in  this 
forecast.   At  the  initial  time  the  two  storms  were  only  660  nm. 
apart.   Studying  a  sample  of  northern  Pacific  storms  over  water, 
Brand  (1970)  found  that  storms  having  separation  distances  of  less 
than  750  nm.  were  subject  to  an  obvious  "Fujiwhara  effect".   Although 
both  storms  had  a  maximum  surface  wind  of  52  knots,  the  maximum 
grid-point  value  of  absolute  vorticity  of  Fern,  216  x  10   sec 

_r  _-i 

was  much  stronger  cyclonically  than  that  of  Edith,  167  x  10   sec 
Edith  being  embedded  in  the  western  periphery  of  the  subtropical 
high  as  noted  earlier.   In  fact,  whereas  Edith  had  become  the  dominant 
vortex  of  the^  two  before  130000Z  (Fern  had  weakened  to  a  tropical  de- 
pression by  120400Z),  the  barotropic  prognosis  maintained  Fern  as 
the  stronger.   Thus  with  Fern  stronger  in  the  forecast,  the  inter- 
action would  have  been  onesided  causing  Edith  to  move  northward  where 
it  was  engulfed  sooner  by  the  strong  westerlies.   Additionally,  the 
mutual  attraction  consistent  with  atmospheric  systems,  but  excluded 
from  theoretical  nondivergent  systems  (Brand,  1970),  could  have  been 
instrumental  in  the  northwestward  track  of  the  real  Edith. 

The  120000Z  forecast  represents  a  re-run.   An  initial  movement 
of  the  forecast  toward  the  southwest  prompted  a  check  of  the  analysis 
winds  nearby.   It  was  immediately  obvious  that  although  Edith  had 
weakened  to  a  35-knot  wind  maximum,  purging  of  the  storm  wind  at 
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Merida,  Mexico,  station  76644  (only  73  ran.  distant)  had  left  a  north- 
westerly recalculated  wind  (335/9.4)  at  Merida.   Therefore,  grid 
points  nearby  naturally  responded  with  northwesterly  winds.   This 
difficulty  was  removed  by  substituting  for  station  76644  a  wind 
which  (when  purged  of  the  model  hurricane  wind)  would  leave  a  flow 
equal  to  the  most  recent  operational  estimate  of  the  storm's  course 
and  speed.   Thus  replacing  the  original  station  76644  wind  108/10.8 
with  125/26  yielded  a  recalculated  wind  of  113/7.7  (close  to  the 
past  12-hour  movement  vector  of  Edith).   The  forecast  using  this 
correction  provided  a  respectable  prognosis  for  only  24  hours  after 
which  the  relative  closeness  of  an  overly  strong  Fern  to  the  west  and 
the  subtropical  high  to  the  east  appeared  again  to  dominate  as  at 
111200Z. 

Forecasts  at  121200Z  and  130000Z  unfortunately  represent  the 
same  effect,  despite  efforts  to  improve  the  121200Z  forecast.   Since 
by  this  time  Fern  had  weakened  to  a  depression,  no  storm  parameters 
were  entered  on  the  original  run  for  that  storm.   Accordingly,  it 
was  theorized  that  since  no  storm  purging  would  be  effected  upon 
the  yet  strong  mean  winds  surrounding  the  vortex  in  southern  Texas, 
the  wind  analysis  for  grid  points  in  the  western  Gulf  would  be  in- 
accurately portrayed.   However  even  with  storm  parameters  entered  for 
Fern,  little  improvement  was  realized. 

Repositioning  the  130000Z  forecast  to  begin  at  the  actual  track 
position,  one  finds  that  the  72-hour  error  of  280  nm.  is  not  excessive 
for  a  sharply  recurving  hurricane. 
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FERN 

Figure  4  displays  the  path  of  Fern  after  it  has  moved  southward 
out  over  the  Gulf.   Also  included  are  the  SANBAR  forecasts  with 
storm  parameters  displayed  in  Table  3.   In  view  of  the  meandering 
path  of  Fern,  the  three  first  forecasts  are  not  poor  ones,  and  the 
090000Z  forecast  comes  within  30  miles  of  the  point  of  landfall  des- 
pite the  90-degree  direction  error. 

Table  3.   Storm  parameters  for  Fern   (Eye 
diameter  =  20  nm.  throughout) 

Forecast  time      Max.  sfc.  wind  (knots)      Influence  radius  (nm.) 

275 
325 
325 
375 
375 
350 
325 
300 

The  error  of  a  northward  forecast  from  091200Z  is  quickly  ex- 
plained when  one  examines  the  recalculated  winds  of  the  surrounding 
stations.   Station  72240  with  an  input  wind  of  140/22.5  had  a  recal- 
culated wind  of  199/11.7.   Likewise,  station  72255  went  from  039/11.6 
to  208/9.9,  and  station  72250  (missing)  had  a  wind  of  227/15.4  com- 
puted from  its  regression  equation.   These  three  stations  naturally 
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Fig.  4  —  Hurricane  Fern  from  0000Z  8  September 
to  OOOOZ  13  September.   For  details  see  Fig.  3. 
(Forecasts  beginning  101200Z  do  not  extend  to 
72  hours.) 
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supplied  a  strong  southerly  component  to  their  surrounding  grid- 
points. 

The  forecast  from  100000Z  was  quite  good  in  predicting  the  point 
of  landfall,  but  again  at  101200Z  an  .immediate  northwestward  movement 
was  predicted.  The  movement,  though  not  completely  understood,  could 
have  been  influenced  by  both  stations  72248  and  72255  having  recalcu- 
lated winds  from  127°. 
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The  shortened  forecasts  from  110000Z  and  111200Z  returned  to  more 
acceptable  errors,  but  by  this  time  Fern  was  moving  inland. 

GINGER 

Figure  5  represents  only  a  small  portion  of  the  record  lifetime  of 
Hurricane  Ginger.   As  before  SANBAR  forecasts  have  been  added  with 
Table  4  listing  the  applicable  storm  parameters. 

Table  4.   Storm  parameters  for  Ginger  (Eye 
diameter  =  20  ran.  throughout) 

Forecast  time      Max.  sfc.  wind  (knots)      Influence  radius  (ran.) 

325 

350 

375 

350* 

350* 

*   indicates  re-run 

The  first  three  forecasts  show  the  typical  slow  bias,  but  the  120000Z 
72-hour  error  is  held  to  180  nm.  when  Ginger  commences  slowing  and 
curving  toward  the  south. 

The  121200Z  forecast  was  a  rerun  for  obvious  reasons.   Due  to  the 
proximity  of  bogus  point  9  (30°N,  60°W)  to  the  hurricane  center  at  the 
time  of  analysis,  12  nm. ,  the  input  wind  for  that  bogus  point  (270/33.0) 
reduced  to  a  recalculated  wind  of  068/17.0.   Accordingly  a  substitute 
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wind  for  bogus  point  9  was  quickly  prepared  from  the  hurricane's 
immediate  past  movement  vector  plus  the  model  hurricane  wind  vector. 
This  new  wind  of  257/58.0  resulted  in  a  recalculated  wind  of  230/10.8 
closely  approximating  the  hurricanes  movement  vector  as  would  have  been 
known  by  the  operational  forecaster.   Full  credit  cannot  be  accepted 
for  the  small  72-hour  forecast  error  of  85  nm.   Not  only  did  a  reverse 
of  the  hurricane's  course  reduce  the  error,  but  the  SANBAR  forecast 
was  forced  toward  the  southeast  by  boundary  limitations  (not  shown). 

The  same  re-run  technique  was  used  on  the  130000Z  forecast  for 
Ginger.   The  substantial  reduction  in  speed  of  this  forecast  can  be 
attributed,  at  least  in  part,  to  the  weakening  of  the  westerly  com- 
ponents at  Ocean  Station  E  as  well  as  at  bogus  points  9  and  10 
(30°N,  50°W). 

HEIDI 

The  forecasts  of  SANBAR  on  tropical  storm  Heidi  (Figure  6  and 
Table  5)  are  presented  mostly  as  a  thing  of  beauty,  noting  that  the 
48-hour  forecast  error  from  130000Z  was  less  than  80  nautical  miles. 
The  forecast  certainly  supports  the  statement  by  Sanders  (1970),  "  .  .  . 
dramatic  acceleration  of  storms  in  the  westerlies  can  be  accounted  for 
baro tropically.". 

Table  5.   Storm  parameters  for  Heidi  (Eye 
diameter  =  20  nm.  throughout) 

Forecast  time      Max.  sfc.  wind  (knots)      Influence  radius  (nm. ) 

121200Z  40  300 

130000Z  45  300 
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Fig.  6  —  Tropical  storm  Heidi  from  1200Z  12  September 
to  OOOOZ  15  September.   For  details  see  Fig.  3.   (Fore- 
casts extend  only  to  48  hours.) 
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DISCUSSION  AND  CONCLUSIONS 

In  reviewing  the  results  of  SANBAR  as  applied  in  the  previous 
section,  it  is  found  that  erroneous  calculations  of  storm  purged 
winds  led  to  poor  forecasts  at  least  twice  for  each  of  the  three 
hurricanes.   The  other  more  nebulous  difficulty  arose  when  two  storms 
came  into  close  proximity.   During  the  infrequent  occurrences  of 
this  latter  difficulty,  it  is  felt  that  operational  forecasters 
must  study  stream  function  and  vorticity  patterns,  and  formulate  sub- 
jective modifications  to  the  forecast  tracks. 

As  for  the  more  prevalent  difficulty  of  calculating  storm  purged 
winds,  two  types  of  corrective  action  are  available.   The  model  storm 
wind  may  be  modified  by  varying  its  input  parameters,  or  the  input 
(observed)  winds  may  be  modified. 

During  this  study,  the  alteration  of  the  storm's  maximum  influ- 
ence radius  (Edith  100000Z)  proved  to  be  advantageous.   Naturally,  the 
application  of  this  technique  may  be  germane  during  operational  situa- 
tions. 

The  several  instances  of  modifying  the  input  wind  observations  near 
the  storm  center  proved  to  be  very  effective  during  this  investigation. 
The  author  is  encouraged  to  have  learned  through  privileged  communica- 
tion that  NHC  will  institute  a  program  modification  during  the  hurricane 
season  of  1972  which  should  practically  eliminate  the  problem  of  the 
recalculated  wind.   This  modified  SANBAR  model,  as  written  by  Dr. 
Arthur  Pike  of  NHC,  requires  that  the  analyzed  wind  at  all  stations 
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within  the  influence  radius  be  represented  by  the  sum  of  the  symmetric 
storm  circulation  vector  and  the  storm  movement  vector. 

In  short,  the  modified  SANBAR  effects  automatically  those  changes 
which  the  author  undertook  at  single  stations  on  only  a  few  forecasts. 
More  importantly,  the  modified  SANBAR  will  apply  these  corrections  for 
every  forecast  and  for  all  stations  within  the  influence  radius.  Com- 
parison of  SANBAR  forecasts  of  this  study  with  those  of  the  modified 
SANBAR  are  displayed  in  Tables  6  and  7. 


Table  6.  Forecast  Error  Averages  (nm.)  for 
Edith  and  Fern  (080000Z,  101200Z, 
110000Z) 

12-hr  24-hr  A 8-hr  72-hr 

SANBAR  (M.I.T.)         60(6)     96(6)    143(6)    161(4) 
SANBAR  (modified)       39(6)     63(6)     95(6)    160(4) 

(Numbers  in  parenthesis  indicate  the  number  of  forecasts 
compared.   All  forecasts  were  repositioned  so  that  their 
initial  positions  agree  with  that  of  the  after-the-fact 
track. ) 


Table  7.   Forecast  Error  Averages  (nm.)  for 
Ginger  and  Heidi  (121200Z  and 
130000Z) 

12-hr  24-hr  48-hr  72-hr 

SANBAR  (M.I.T.)         44(4)     80(4)    171(4)     95(2) 
SANBAR  (Modified)       29(4)     59(4)    188(4)    272(2) 

(Numbers  in  parenthesis  indicate  the  number  of  forecasts 
compared.   All  forecasts  were  repositioned  so  that  their 
initial  positions  agree  with  that  of  the  after-the-fact 
track.) 
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In  defense  of  this  small  and  selective  sample,  it  must  be  stated 
that  results  for  the  modified  SANBAR  were  available  only  for  the  5 
times  displayed.   Recognizing  that  the  sample  size  is  inadequate  as 
well  as  that  input  storm  parameters  and  bogus  winds  differed  for  the 
two  forecast  systems,  only  during  Ginger  and  Heidi  did  the  original 
SANBAR  have  smaller  errors  than  did  the  modified  SANBAR.   In  partic- 
ular during  both  original  SANBAR  forecasts  for  Ginger,  corrections 
were  made  for  recalculated  winds.   Indirectly, too,  this  may  indicate 
the  worth  of  using  satellite  film  strips  in  preparing  oceanic  bogus 
data,  for  it  is  here  that  the  original  SANBAR  did  better  at  48  and 
72  hours. 

In  regard  to  the  slow  bias  of  SANBAR  which  has  become  well  known 
to  the  operational  forecasters,  subjective  modification  must  be 
attempted  awaiting  such  time  as  an  objective  improvement  can  be 
offered.   A  suspected  equivalent  barotropic  effect  (NHC,  1972) 
which  tentatively  suggests  multiplying  the  individual  winds  by  about 
1.3  has  not  been  successful  to  date. 

In  addition  to  a  possible  new  baroclinic  model  (NHC,  1971),  im- 
provements to  the  statistical  and  analog  methods  as  well  as  those  to 
SANBAR  should  foretell  of  a  more  promising  future  for  the  hurricane 
track  forecaster. 
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